Temperature dependence of the paramagnetic spin excitations in BaFe 2 As 2 
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We use inelastic neutron scattering to study temperature dependence of the paramagnetic spin ex- 
citations in iron pnictide BaFe2As2 throughout the Brillouin zone. In contrast to a conventional local 
moment Heisenberg system, where paramagnetic spin excitations are expected to have a Lorentzian 
function centered at zero energy transfer, the high-energy (hcj > 100 meV) paramagnetic spin exci- 
tations in BaFe2As2 exhibit spin- wave-like features up to at least 290 K (T = 2.1T/v). Furthermore, 
we find that the sizes of the fluctuating magnetic moments (m 2 ^ « 3.6 fi% per Fe are essentially 
temperature independent from the AF ordered state at 0.05T/V to 2.1T/v, which differs considerably 
from the temperature dependent fluctuating moment observed in the iron chalcogenide Fei.iTe [I. 
A. Zaliznyak et a/., Phys. Rev. Lett. 107, 216403 (2011).]. These results suggest unconventional 
magnetism and strong electron correlation effects in BaFe2As2. 



PACS numbers: 75.30.Ds, 25.40.Fq, 75.50.Ee 

The elementary magnetic excitations (spin waves and 
paramagnetic spin excitations) in a ferromagnet or an an- 
tiferromagnet can provide direct information about the 
itinerancy of the unpaired electrons contributing to the 
ordered moment. In a local moment system, spin waves 
are usually well-defined throughout the Brillouin zone 
and can be accurately described by a Heisenberg Hamil- 
tonian in the magnetically ordered state. The total mo- 
ment sum rule requires that the dynamical structure fac- 
tor S(q,u), when integrated over all wave vectors (q) 
and energies (E = TiuS), is a temperature independent 
constant and equals to (m 2 ) = (g/i#) 2 S{S + 1), where 
g is the Lande g factor (« 2) and S is the spin of the 
system pQ. Upon increasing temperature to the para- 
magnetic state, spin excitations in the low-g limit can 
be described by a simple Lorentzian scattering function 
S(q,u>) oc [l/(«i +g 2 )][r/(r + cj 2 )], where n x is the tem- 
perature dependent inverse spin-spin correlation length 
and T is the wave vector dependent characteristic energy 
scale [1H4]. At sufficiently high temperatures above the 
magnetic order, spin excitations should be purely para- 
magnetic with no spin- wave- like correlations. Therefore, 
a careful investigation of the wave vector and energy de- 
pendence of spin excitations across the magnetic ordering 
temperature can provide important information concern- 
ing the nature of the magnetic order and spin-spin cor- 
relations. For example, a recent inelastic neutron scat- 
tering study of spin excitations in one of the parent com- 
pounds of iron-based superconductors, the iron chalco- 
genide Fei.iTe which has a bicollinear antiferromagnetic 
(AF) structure and Neel temperature of Tn = 67 K [5]- 
[TT] , reveals that the effective spin per Fe changes from 
S ~ 1 in the AF state to S ~ 3/2 in the paramagnetic 
state, thus providing evidence that Fei.iTe is not a con- 



ventional Heisenberg antiferromagnet but a nontrivial lo- 
cal moment system coupled with itinerant electrons [12] . 

Since antiferromagnetism may be responsible for elec- 
tron pairing and superconductivity in iron-based super- 
conductors [13l [14] , it is important to determine if the 
observed anomalous spin excitation behavior in iron tel- 
luride Fei.iTe is a general phenomenon in the parent 
compounds of iron-based superconductors. For this pur- 
pose, we study the spin excitations of another parent 
compound of iron-based superconductors, the iron pnic- 
tide BaFe2As2 which has a collinear AF structure with 
T/v ~ 138 K [15-18 , over a wide temperature range 
(0.05T/V < T < 2.1T/v). In the low-temperature or- 
thorhombic phase (T = 0.05T/V = 7 K), previous inelastic 
neutron scattering experiments found that a Heisenberg 
Hamiltonian with highly anisotropic effective magnetic 
exchange couplings and damping along the orthorhom- 
bic a and b axes directions can describe the observed 
spin- wave spectra [19 , similar to the spin waves in the 
collinear AF ordered CaFe2As2 [20]. However, similar 
measurements on iron pnictide SrFe2As2 suggest that 
spin waves can be better described by calculations from a 
five-band itinerant mean- field model [2TJ[22]. Therefore, 
it is unclear whether a localized Heisenberg Hamiltonian 
[T9] [20] or itinerant magnetism [21- 23 is a more appro- 
priate description for spin waves in pnictides. In the high- 
temperature tetragonal phase, the spin-wave anisotropy 
of BaFe2As2 appears to persist at 150 K (T = l.lT/v) 
suggesting the presence of an electronic nematic phase 
[T9] [241129] . However, these paramagnetic spin excita- 
tions can also be understood by considering both the lo- 
calized and itinerant electrons using dynamic mean field 
theory [30] or a biquadratic spin-spin interactions within 
a Heisenberg Hamiltonian without the need for electronic 
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nematicity [3T| [32]. 

By studying spin excitations in BaFe2As2 over a wide 
temperature range throughout the Brillouin zone in abso- 
lute units, we can determine the temperature dependence 
of the paramagnetic scattering and its spectral weight. 
This will reveal if itinerant electrons in BaFe2As2 are 
coupled with local moments on warming across TV sim- 
ilar to that of the iron telluride Fei.iTe [12 . Surpris- 
ingly, we find that the total fluctuating magnetic mo- 
ments (m 2 ) ~ 3.6 n 2 B per Fe in BaFe2As2, correspond- 
ing to an effective spin 5 = 1/2 per Fe [33], are essen- 
tially unchanged on warming from 7 K at T = 0.05T/V to 
room temperature at 2.1T/v, much different from that of 
Fei.iTe [12]. In addition, while paramagnetic spin exci- 
tations at small wave vectors near the AF zone center fol- 
low a simple Lorentzian scattering function as expected 
[2], they change only slightly from the low-temperature 
spin waves for wave vectors near the zone boundary up 
to room temperature. This is different from the expecta- 
tion of a local moment Heisenberg system, and indicate 
a strong electron correlation effect in BaFe2As2- 

We have used the MAPS time-of-flight inelastic neu- 
tron spectrometer at ISIS, Rutherford- Appleton Labora- 
tory, UK, to determine the paramagnetic excitations of 
BaFe2As2- For the experiment, we have used the same 
sample and experimental set-up as described previously 
[19]. Below T/v, BaFe2As2 has an orthorhombic struc- 
ture with a = 5.62 and b = 5.57 A and forms a collinear 
AF order at the ordering wave vector Q = (1,0, 1) [18] . 
In the paramagnetic state, BaFe2As2 changes to tetrag- 
onal structure. Figure 1 presents an overview of the 
temperature evolution of the spin excitations at differ- 
ent energies. The data has been normalized to a vana- 
dium standard and plotted in absolute units of mbarn 
gr -i me V _1 f.u. -1 , without correction for the magnetic 
form factor, leading to a decrease in magnetic scatter- 
ing with increased Q. At E = 50 ± 10 meV, spin 
waves form ellipses along the transverse direction cen- 
tered at Qaf = (1 + n, L) and Qaf = 1 + ^ L), 
where m, n = 0, 1, 2, • • • and L = 1, 3, 5, • • •, at 7 K (Fig. 
la) [19] . Upon warming to the paramagnetic state at 
225 K (T = 1.63TV, Fig. lb) and 290 K (2.17V, Fig. 
lc), the signal becomes weaker, and the ellipses broader, 
compared to the spin wave peak seen at 7K, similar 
to the low-energy paramagnetic spin excitations seen in 
CaFe2As2 [34]. However, the spin waves at E = 100 ± 10 
(Figs, ld-lf) and 150 ± 10 meV (Figs, lg-li) only de- 
crease slightly in intensity on warming, and become more 
diffusive at 290 K. 

Figures 2a-2d show the background subtracted scat- 
tering for the Ei = 450 meV data projected in the wave 
vector (Q = [l,if]) and energy space at T = 7, 125, 
225, and 290 K, respectively. The solid lines represent 
the expected dispersion from the anisotropic Heisenberg 
Hamiltonian [19 . At 7 K, three plumes of spin waves 
stem from Qaf = [1, K] where K = 0, ±2, and reach to 
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FIG. 1: (Color online) Panels a)-c) compare the E — 50meV 
magnetic scattering deep inside the ordered state (7 K) to 
scattering in the paramagnetic phase for temperatures (225 K 
and 290 K) well away from the TV = 138 K phase transition. 
Panels d)-f) and g)-i) are a similar comparison for energy 
transfers of 100 meV and 150 meV, respectively. The dotted 
ellipses and boxes are guides to the eye to more easily facilitate 
comparison. Data in panels a)-f) and g)-i) were collected 
using Ei — 250 meV and 450 meV respectively. All data 
was background subtracted using the average intensity from 
the region 1.8 < H < 2.2, -0.2 < K < 0.2 r.l.u. as the 
background point. Data in the region H < was folded into 
the equivalent H > positions in order to improve statistics. 



the zone boundary at ~200 meV (Fig. 2a). On warming 
to 125 (Fig. 2b), 225 (Fig. 2c), and 290 K (Fig. 2d), spin 
excitations become broader in momentum space but their 
zone boundary energies appear to be unchanged. For 
the classical insulating Heisenberg ferromagnet or anti- 
ferromagnet, spin excitations in the paramagnetic state 
should be uncorrelated and display Lorentzian-like peaks 
centered at E = meV at sufficiently high temperatures 
[2] [3]. If electron correlations are important, spin excita- 
tions in the paramagnetic state should exhibit spin-wave 
like peaks in energy for wave vectors near the zone bound- 
ary [4] . While previous work found that spin excitations 
near the zone boundary for energies above E = 100 meV 
are indeed similar between 7 K and 150 K [19 , it is un- 
clear what happens to zone boundary spin excitations at 
higher temperatures. 

Figure 3 summarizes the wave vector and temperature 
dependence of the spin excitations from 7 K to 290 K 
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FIG. 2: (Color online) Panels a)-d) compare the Energy 
versus K intensity slices for 7 K, 125 K, 225 K, and 290 K. 
All data is background subtracted and folded in an identical 
manner as described in the caption of Fig. [I] The solid 
line is the Heisenberg dispersion obtained using anisotropic 
exchange couplings SJ la = 59.2 ± 2.0, SJ lb = -9.2 ± 1.2, 
SJ 2 = 13.6 ± 1.0, SJ C = 1.8 =b 0.3 meV determined by fitting 
the full cross section to the 7 K data [19]. 
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FIG. 3: (Color online) a-f) Temperature overplot of the evolu- 
tion of the spin excitations as a function of increasing energy. 
The green diamonds, yellow squares, red circles, cyan upward 
facing triangles, and blue downward facing triangles are for 
the 7 K, 125 K, 150 K, 225 K, and 290 K data respectively. 
The data has been artificially offset for clarity and empirically 
fit using Gaussian functions. The insets are the fits without 
offset. 



along the Q = [1,K, 0] direction. For each of the wave 
vector cuts along the if-direction, the if-direction inte- 
gration range is slightly different. At 10 < E < 20 meV, 
the spin wave intensity increases on warming from 7 K to 
125 K. Upon further warming to above T/v, the spin exci- 



FIG. 4: (Color online) a) Dispersion along the [1,K] direc- 
tion as determined by energy and Q cuts of the raw data. 
The solid line is the anisotropic Heisenberg dispersion [19] . 

b) Dispersion along the [if, 0] direction built using the same 
method. The light blue upward facing triangular points in 

c) -f) are constant-Q cuts at Q = (1,0.05), (1, 0.2), (1, 0.35), 
and (1, 0.5), respectively, at 225 K. The dark blue downward 
facing triangular points in c)-f) are identical constant-Q cuts 
at 290 K. The solid green and red lines are guides to the eye 
describing the observed 7 K and 150 K scattering, respec- 
tively. These constant-Q cuts correspond to cuts across the 
dispersion as depicted in the inset of panel c) . The horizontal 
bars in d) and f) are instrumental energy resolution. 

tation peak centered at Qaf = (1,0, L) becomes weaker 
and broader with increasing temperature, and is very 
broad at 290 K. For spin wave energies 20 < E < 30 meV 
(Fig. 3b), 30 < E < 40 meV (Fig. 3c), and 40 < E < 50 
meV (Fig. 3d), the situation is similar although spin exci- 
tations have less temperature dependence with increasing 
energy. Finally, spin excitations only change marginally 
from 9 K to 290 K for 90 < E < 100 meV (Fig. 3e) and 
125 < E < 145 meV (Fig. 3f). 

Based on the data in Figures 2 and 3, we construct in 
Fig. 4a spin excitation dispersions along the [1,K] and 
[if, 0] directions at 225 K (the upper triangles) and 290 
K (the lower triangles). Comparing the outcome with 
the spin waves at 7 K (the solid lines) reveals essentially 
the same dispersion for spin excitations in the param- 
agnetic state for temperatures up to T = 2.1T/v (Figs. 
4a and 4b). Figures 4c-4f show constant-Q cuts of the 
spin excitations along the [1, K] direction throughout the 
Brillouin zone (see inset in Fig. 4c). Previous measure- 
ments at 7 K and 150 K are plotted as green and red solid 
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FIG. 5: (Color online) The local susceptibility plots in pan- 
els a)-e) represent the total Q-integrated intensity across the 
magnetic zone of size < H < 2, — 1 < K < 1 r.l.u. In prac- 
tice, it is not possible to use the actual full zone size in H and 
K because of gaps in the detector array and consequent lim- 
ited accessibility of certain reciprocal space regions. Thus, for 
each Ei a smaller region that either contains all of the scat- 
tering and/or has the requisite symmetry is chosen instead. 
These regions are then all normalized to the entire zone area 

J X "(Q,E)dQ 

in panels a)-e) are empirical fits of the local susceptibility, f) 
An overplot of these fits to aid in a cross-comparison of the 
temperature dependence, g) The dynamic moment as deter- 
mined by integrating the fits from the previous panels. The 
static moment is reproduced from [18] . 



lines, respectively. For wave vector near the zone center 
at Q = (1, 0.05) and (1, 0.2), we see that the well-defined 
spin wave peaks in the AF phase become Lorentzian like 
in the paramagnetic state at 150 K (Fig. 4c and 4d). On 
further warming to 225 and 290 K, quasielastic intensity 
near E — meV becomes weaker, consistent with the ex- 
pectations for paramagnetic scattering [2H1]. However, 
the low-temperature spin wave peaks at the wave vec- 
tors Q = (1,0.35) and (1,0.5) near the zone boundary 
are still clearly present up to 290 K, and only become 
slightly broader and weaker (Figs. 4e and 4f), thus sug- 
gesting a strong electron correlation effect in BaFe2As2- 



Finally, we show in Figure 5 the temperature depen- 
dence of the local dynamic susceptibility for BaFe2As2 
[33] [35]. In the AF ordered state at 7 K, there is a spin 
anisotropy gap below ~10 meV [36] and the local sus- 
ceptibility peaks at ~180 meV (Fig. 5a). On warming 
to 125 K just below T/v, the spin anisotropy gap dis- 
appears while at higher energies the local susceptibility 
remains essentially unchanged (Fig. 5c). Upon further 
warming to the paramagnetic state at 150 K (Fig. 5e), 
225 K (Fig. 5b), and 290 K (Fig. 5d), we see that the 
local dynamic susceptibility becomes slightly weaker and 
broader with increasing temperature (Fig. 5f). Figure 
5g shows the temperature dependence of the ordered mo- 
ment (solid line) [18] and integrated local susceptibility, 
which is dominated by spectral weight from spin exci- 
tations above 100 meV. For comparison, we note that 
the integrated magnetic spectral weight of Fei.iTe were 
reported to concentrate almost entirely within 30 meV 

In earlier triple-axis spectrometry studies of paramag- 
netic spin excitations of metallic ferromagnets such as 
iron and nickel, there was considerable controversy con- 
cerning whether persistent spin wave like excitations can 
exist in the paramagnetic state above Tq [3] I37H40] . For 
BaFe2Aa2, we see spin- wave- like excitations above 100 
meV at temperatures up to 2.1T/v- This is different from 
the usual paramagnetic scattering in a Heisenberg anti- 
ferromagnet. The lack of temperature dependence of the 
integrated local moment, (m 2 ) « 3.6 fi B per Fe, suggests 
that the effective spin of iron in BaFe2As2 (S — 1/2) 
is unchanged from the AF orthorhombic phase to the 
paramagnetic tetragonal phase up to room temperature. 
Therefore, there is no exotic entanglement of itinerant 
electrons with localized magnetic moments, much differ- 
ent from that of the Fei.iTe [12]. We also note that the 
size of the dynamic moment (m 2 ) « 1.9 \±b per Fe in 
BaFe2As2 is larger than the local moment of 1.3 \±b per 
Fe determined from x-ray emission spectroscopy [41] , but 
similar to the local moment of 2.1 \ib per Fe in SrFe2As2 
obtained from the Fe 3s core level photoemission spectra 
measurements [42] . 

In summary, we have studied the temperature de- 
pendent paramagnetic spin excitations in iron pnictide 
BaFe2As2, one of the parent compounds of iron-based 
superconductors. In contrast to a conventional Heisen- 
berg system, we find spin-wave-like paramagnetic exci- 
tations near the zone boundary for temperatures up to 
2.1T/v with no evidence for the expected zone bound- 
ary magnon softening. In addition, the integrated local 
magnetic moment is remarkably temperature indepen- 
dent from the AF ordered orthorhombic phase to the 
paramagnetic tetragonal phase, and corresponds to an 
effective iron spin of S = 1/2. This is different from 
the temperature dependent spin excitations in the iron 
chalcogenide Fei.iTe. Our results indicate a strong elec- 
tron correlation effect and suggest that the entanglement 
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of itinerant electrons with localized magnetic moments in 
Fei.iTe [12] is not fundamental to the magnetism in the 
parent compounds of iron-based superconductors. Fur- 
thermore, correctly modeling the pnictides requires tak- 
ing into account a mixed state where correlations are 
important. Indeed, both dynamic mean field theory [30] 
and biquadratic exchange [3TJ [32] are approaches that 
pick up electron correlations and appear to provide neces- 
sary features for describing the physics of these systems. 
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